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Abstract 

This article gives a comprehensive review on the recent studies of the double exchange 
systems using non-perturbative approaches; the dynamical mean-held theory and the Monte 
Carlo method. Investigations beyond mean-field type treatments are described. Taking into 
account strong spin fluctuations which create large changes in conduction electron structure, 
finite temperature properties as well as dynamics of the system are calculated. Comparisons 
with experimental data for colossal magnetoresistance manganites are made. We show that 
high Curie temperature (T c ) compounds, e.g. (La,Sr)Mn03, are canonical double-exchange 
systems. Properties of other compounds with lower T c are discussed in relation to inhomo- 
geneities of the system including the issue of phase separation. 



I. INTRODUCTION 

The concept of the double-exchange (DE) interaction was introduced by Zener& in order to 
explain the ferromagnetism of the perovskite manganites AM.Y1O3. He considered a Kondo-lattice 
type model 

H = ~ { C ia C jv + h - C ) ~ Ju^2^i- Si, (1) 

ij,a i 

where t and Jh > are e g electron's hopping and (ferromagnetic) Hund's coupling between e g 
and ^2g electrons, respectively. An effective Hamiltonian in the limit Jh — ► 00 was introduced by 
Anderson and Hasegawa& in the form 

n = -Y J t{s l1 s 3 ){c\~c J + h.c). (2) 

Mean-field type arguments of these models, including those by de Gennes,!) helped us to discuss 
the magnetism of manganite compounds. 

However, recent reinvestigations of these models which are motivated by the observation of 



colossal magnetoresistance (CMR) phenomenaa EP in manganites revealed that such simplified 
treatments are insufficient to discuss quantitative nature of the models. Let us show an example. 
Curie temperature T c of the model has been estimated by Millis et al& using the mean-field 
type discussion. Using some appropriate values for t and Jh, they found that the mean-field T c 
of the model is much larger than those for perovskite manganites (La,Sr)Mn03, and concluded 
that double-exchange alone is insufficient to explain the thermodynamics of these manganites. 
However, as we will show in this article, accurate treatments for the model give suppression of T c 
from the mean-field value, and show consistency with (La,Sr)Mn03 data. 
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Another important point which will be discussed here is the systematic understanding of ex- 
periments. For example, it is known that CMR manganites exhibit metal-insulator transition at 
Curie temperature. This statement is, however, inaccurate. By varying compositions, varieties of 
phases with different properties have been knownJZP A typical reference compound La2/3Sri/3Mn03 
shows no metal-insulator transition. The resistivity always increases by the increase of tempera- 
ture, i.e. dp(T)/dT > even above T c . Its absolute value in the paramagnetic phase is around the 
Mott's limit value, and may be explained by the DE model. On the other hand, La2/3Cai/3Mn03 
which is often reffered to as an "optimal" CMR material shows insulating behavior above T c with 
much larger p(T). Thus we have to be mindful to make distinctions between various compositions 
of manganites in a systematic manner. 

The purpose of this article is to solve such misunderstandings and confusions in both the- 
oretical MJid ,£Xfte^r^ntal studies. We review the investigation of the DE model by the au- 
thor.&&ll3llJ'0'll3llJ'y'll3* The first part of this article shows the finite temperature behaviors in 
the DE systems. In §2 we introduce the model and the method (dynamical mean-field approach 
and the Monte Carlo calculation). In §3 the results for infinitely large spin (classical spin limit) is 
presented, and the 1/5* correction is shown in §4. The second part is devoted for discussions with 
respect to comparisons of model behavior with experimental data. Comparison with experiments, 
are shown in §5. Section 6 is devoted for concluding remarks. 



II. MODEL AND METHODS 

2.1 Double Exchange Model in the Large Spin Limit 

Model. The compound LaMn03 has four 3d electrons per atom in the (t2g) 3 (e g ) config- 
uration. Due to Hund's coupling, these electrons have the high spin state, i.e. spin parallel 
configuration. By substitution of the La site with alkaline-earth divalent ions, holes are doped as 
carriers which is considered to enter the e g orbitals. 

As briefly mentioned in the previous section, Zener& introduced a Kondo-lattice type Hamilto- 
nian ([!]) with ferromagnetic spin exchange Jh between localized spins and itinerant electrons. In 
manganites (R,A)Mn0 3 , Hund's coupling Jh is estimated to be a few eV while electron hopping t 
be in the order of O.leV. Therefore, we have to deal with a strong coupling region Jh ^ t. In order 
to discuss the ferromagnetism of the model in the metallic region, Zener introduced the notion of 
the "double exchange" interaction. 

In some cases, the limit Jh — ► oo first studied by Anderson and HasegawaEP is called the DE 
model. However, in order to avoid complications, we call the model (|XJ) in the strong coupling 
region as DE model. In the weak coupling limit Jn/t <C 1, the model is often referred to as the 
s-d model and studied as a model hamiltonian for the magnetic semiconductors. 

The classical rotator limit, or equivalently large spin limit S = oo, has been introduced by 
Anderson and Hasegawa.0) In manganites, we consider the case where the localized spin is in a 
high-spin state (S = 3/2) with the ferromagnetic coupling, the effect of quantum exchange seems 
to be less relevant compared to thermal fluctuations, at least in low energy physics. However, 
the role of quantum exchanges might give non-trivial effects in the system in the region with less 
thermal fluctuations. Such issues are left for future studies. 

In this article, we introduce a model with finite Jh and infinite S, 



Hereafter we express the localized (classical) spin by m; = (rrii x ,rrii y ,rrii z ) with the normalization 




(3) 



rn 
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Previous investigations. Anderson-Hasegawall as well as de Gennes& studied the mag- 
netism of the model with localized spins treated as static, i.e. neglecting spin fluctuations, they 
made a search for energetically favored spin configurations as well as a mean-field calculation for 
finite temperatures. However, if we consider the system at finite temperature, especially around 
T ~ T c , spin fluctuations 8 Si become important. Especially, in our case of Jh 3> t, such spin 
fluctuations give large effects to electronic structures since J^SS 3> t. 

The point of interest for us in relationship with CMR phenomena is the change of resistivity 
p at around T c . Mean-field type treatments axe not justified for this purpose. An alternative 
approach has been made by Kubo and Ohata.E-r They phenomenologically assumed the electron 
self-energy £ in the form 

r _1 = 1m £ oc 1 — M 2 , (4) 

where r is the quasiparticle lifetime, and calculated the resistivity using the Drude's formula 
p oc r -1 . Within this phenomenological treatment, the result qualitatively helps us to understand 
the magnetoresistance (MR) of the DE model via spin disorder scattering mechanism. However, 
in a quantitative way, it fails to reproduce the MR of (La,Sr)Mn03 as we will show in §iTI. 



Beyond the previous theories. Thus, in order to understand the behaviour of the DE 
model and to make direct comparisons with experimental data, it is very important to treat the 
model in an accurate way. The methods have to be an unbiased non-perturbative approaches since 
Jh ^> t, and must take into account the effect of spin fluctuations to calculate finite temperature 
behaviors. 

Here we introduce two methods: The dynamical mean-field (DFM) theory and the Monte Carlo 
calculation. Both of these methods give us the electronic states at finite temperature including 
T ~ T c . They are unbiased and become exact in the limit of large coordination number and large 
system size, respectively. One of the advantages of these methods, in the viewpoint of comparison 
with experiments, is that it is easy to obtain dynamical quantities such as density of states (DOS) 
and optical conductivity a{uj). 

2.2 Dynamical mean-field theory 

We introduce the DMF method for the double-exchange model. For a general review of the 



field, see the review articles in refs. [18] and |19| Within the general scheme of the DMF, we treat 
a lattice system by considering a single-site coupled with "electron bath" , or the time-dependent 
mean-field Gq. This method becomes exact in the large coordination number limit or equivalently 
large spatial dimension limit .E2r ElP 

Generic part of the DMF treatment is as follows. Solving a model-specific single-site problem, 
we obtain the self-energy T<(iLd n ) from Gq. Lattice Green's function is approximated by 

G(k, w n ) = [iw n -{e-y)- S(iu;„,)] _1 . (5) 
The local Green's function is defined by 

Gloc(i^n) = J7 G ( k ' ( 6 ) 
k 

Since X is /^-independent, k dependence of G(k, iu n ) comes through the energy dispersion and we 
have 

G loc (iu n ) = f deN {e)[iu n - (e - n) - t{iu n )]-\ (7) 



N (e) is the DOS for the noninteracting lattice system. The information of the lattice geometry 
is included through the noninteracting DOS in eq. (0). The method is applicable to finite size 
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systems by taking the sum over k points in eq. (H) in the discrete /c-space. We self-consistently 
obtain the time-dependent mean-field G as 

G (iuj n ) = (Gy l(iuj n ) + t(iu n ^ 



Let us now discuss the model-specific part.1) For the present system the action of the 
effective single-site model is described as 



S = -f dnf dra^CrOGo^ri-ra)*^) 
Jo Jo 

-J H / drm-r T fff t . (9) 
Jo 



Here = (cJ, cj) is the Grassmann variables in the spinor notation. Green's function in the 
imaginary time is calculated as 



G{iu n ) = (^(G 1 (iuj n ) + J u rh ■ aj 

= I [ dfiP(m) (Gq 1 ^) + J H m -aY 1 . (10) 

Zj J 

Pirn) is the Boltzmann factor for the spin 

P(m) = lexp[-S cS (m)}, (11) 

Zj 

where S e s is the effective action for the spin 

S e g(m) = -logTr F exp(-S) 

1 e iuJn0+ . (12) 



5^ log det ^—{Gq 1 (iu n ) + J n rha) 



Z is the partition function 

Z = J dn^ J V^*Vty exp(-S) 

= J d^exp[-S eff (m)]. (13) 

The self energy for the single-site system £ is obtained from E(io; n ) = Gq (iu n ) — G _1 (ic<j n ). 
Magnetization of the local spin is obtained by 

(m) = / dQ m P{m)fn. (14) 



Hereafter we take the axis of the magnetization in z direction and the order parameter is expressed 
as M = (m z ). Transport .properties are obtained through the Kubo formula. Conductivity in 
D = oo is calculated aBBEi 

criuj) = (To / I duj IAuj ,u> + uj) , (15) 



where 



Ulu u lu 2 ) = I N (e)de W^A^e, LO X )A a {e, u 2 ). (16) 



Here, the spectral weight function is defined by 



1 

A tr (e,u) = — lmG a (e,uj + irf), (17) 

7T 

while / is the Fermi distribution function. The constant <jq gives the unit of conductivity. In 
this formula, we used that the vertex correction cancels in the conductivity calculation at in 
dimensional limit.^ Thermopower is also calculated: The Seebeck coefficient is obtained as 



28, 18D 



where (k — 1, 2) is defined by 



L h = E/^ (-^) i^oOG&O*- 1 - (19) 

The method is easily expanded to a Bethe lattice with two-sublattice symmetry. In this case, 
magnetic phases with ferromagnetic and antiferromagnetic order parameters can be considered 
simultaneously. We introduce a = A, B sublattice indices for the Weiss fields Go a (iuj n ), and solve 
coupled self-consistency equations. The formula to calculate the Green's function is now given by 

G a {iuj n ) = J dQ a P a {m) (G^(iu n ) + Jm-ay\ (20) 

Here, the Boltzmann weight for the configuration of local spin P a (rn) is calculated from the 
effective action S a , 

§ a (m) = -logTr F exp[-S(G 0cn m)], (21) 

P a {m) = exp[-S a {m)}/Z a , (22) 

Z a = dO a exp[— S a (fh)]. (23) 



Integration over DOS as in eq. @ gives the self-consistent mapping relationil) 

G A _1 (iw n ) = \uj n + n-G B (iun) W 2 /4, 

G 0B -\iu n ) = iu n + fi-G A (iu n )W 2 /4. (24) 



Now the self-consistency equations (|20|) -(24) form a closed set. Within this approach, we can study 
the instability and the formation of magnetic ordering with ferromagnetic, antiferromagnetic, and 
canted antiferromagnetic symmetries. 

2.3 Analytical solutions of the dynamical mean- field theory 

In some limiting cases, analytical solutions are available for the DMF calculations. Analytical 
expressions are in general quite useful to obtain intuitions. 

Paramagnetic phase. For the paramagnetic solution, Green's function is a scalar function 
with respect to spin rotation, so we have G (iuj n ) = g (iu n )I where I is the 2x2 eigenmatrix. 
Then, we have 

G(i^) = ( § °' i "i? "/"% ). (25) 
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Since (m) = and (|^| 2 ) = 1, we have 

G(iu n ) = — — — — I = - — — + — — I. (26) 

# (l^n) - Ju 2 2 \g (lU n ) 1 + J H Po(l^n) ^Jh/ 

The self-energy is then given by 

£(iuv) = G?o - G~ l {\u) n ) = J n 2 G {\u n ). (27) 

From the derivation, we see that the Green's function (|26|) is the same as that of the system 
with the Ising substrate spin in = (0,0, ±1). Furthermore, the Green's function (pBj) shares the 
same analytical structure as that of the infinite-dimensional Falicov-Kimball model (FKM) (or, 
simplified Hubbard model) 

n FK M = - E tiM c i + h - c -) + U Y, c tafifi (28) 



at (rif) = 1/2. Green's function of the FKM in infinite dimension is described ascS'™ 

G FKM (i^) = + . (29) 



Then, the Green's functions ( |26| ) and ( |29| ) share the same analytical structure at {rif) = 1/2. 
In the FKM, the c-electrons are scattered by the charge fluctuations of the localized f-electrons, 
which corresponds to the scattering process of the itinerant electrons by the localized spins in the 
DE model. 

Hence, thermodynamical properties of the DE model in D = oo and S = oo can be understood 
from the nature of the FKM in D = oo which has been studied intensively. For example, the 
imaginary part of the self-energy at the fermi level is finite XmS(O) 7^ in the paramagnetic 
phase.il'ii' 

In Fig. |I|(a) we schematically illustrate the spectral function A(k, uj) in the paramagnetic phase. 
As in the case for FKM, the spectral weight is split into two parts at uj ~ ± Jh for sufficiently large 
Hund's coupling Jh ^ W. For_the semi-circular density of states with the bandwidth W, we have 
the metal-insulator transitionEf at Jh° — 0.5PF, whicLhas the Hubbard-III like nature. Kondo 
resonance peak, which is seen in the Hubbard modeled SB* is missing in this model, since the 
quantum exchange process is absent. As the magnetic moment is induced, the imaginary part of 
the self-energy decreases because the thermal fluctuation of spins decreases. At the ground state, 
spins are magnetically ordered. For the ferromagnetic ground state, there exists two free electron 
bands which are energetically split by Jh exchange interactions. 

Infinite Jh limit —Green's function—. In the case of Lorentzian DOS, Green's function 
is easily obtained in the limit Jh — > 00 Jl^P We consider the hole doped region where \i ~ — Jh- In 
the case of the Lorentzian DOS, the self-consistency equation gives 

G (u + ir]) = (VI- J R + iW)-\ (30) 

Here, chemical potential is /i = — Jh + Sfi where 5fi = 0(W), and Q = to + 5fi = 0(W) is the 
energy which is measured from the center of the lower sub-band — Jh- 

Magnetic field in the z direction is applied to the localized spins in the paramagnetic phase, 
and the induced magnetization is expressed as M = (m z ). Since G in eq. (|30| ) is spin independent 
even in the spin polarized cases, eqs. (|25| ) and (^) gives 

G a (oj + 177) = — — . rT - 

(Q- J H + iW) 2 - J H 2 

1 + Ma 1 



Q + iW 



+ 0(1/J H ). (31) 
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M=0 



0<M<1 



M=1 






(a) 



(b) 



(c) 



Figure 1: Schematic behavior in the spectral function for up-spin electrons, at (a) param- 
agnetic state M = 0, (b) ferromagnetic state at < T < T c , and (c) at the ground state 
T = and M = 1. Solid curves illustrate peak positions of A(k,uj). Width of the curves 
represent height of the peak (quasiparticle weight), while error bars represent the linewidth 
(inverse of lifetime). Grey lines are guides to eyes. Spectral functions for down-spin electrons 
are obtained by exchanging upper part (u ~ Jh) and lower part (uj ~ — Jh). Also see Fig. ||| 
in SI3.1I for actual data. 



At Jh/W — > oo, the spectral weight is calculated as 

* / \ 1 _ , . , 1 + Ma 1 W 

M») = --ZnG.{u, + v 1 ) = —j-^ lWTm . (32) 

We see that the center of the spectral weight is indeed shifted to — Jh- The amplitude of A a is 
proportional to the population of the local spins parallel to a, which indicates that the electronic 
states that are anti-parallel to the local spin are projected out. The self-energy is calculated from 
eqs. © and (13) as 

E^u + irj) = _J H -i^jg(n + iW). (33) 



Eq. (j33|) gives 7?eS ~ — Jh, so the shift in [i is self-consistently justified again. 

Similarly, Green's function at uj' ~ 2Jh, namely at the upper subband, is described as follows. 
Using uj' = uj + 2Jh, where uj = 0(W), Green's function is given by 

G a (u + 2J H + irj) = l ~z\l. lW + 0(1/ J*), (34) 
and the spectral weight is calculated as 



From eqs. ([32]) and (]35), we see the transfer of the spectral weight by magnetization. 
From above equations, Green's function is given in the form 



G a (k,u; M) 



z«(M) 



uj+J R + fi- d l) (k; M) + irS } (M) 

z^(M) 

uj-Jn + fi- (i u \k; M) + ir^(M) 



+ ~ *M,u. ^ i ;rWr- ' (36) 
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asymptotically at Jn/W — > oo. In this limit, Green's function is a sum of contributions from 
lower subbandrs (1) and upper subband (u). Quasiparticle residue is given by 

z®{M) = P CT + (M), 4 U )(M) = P-(M), (37) 

quasiparticle dispersion relation is described as 

C®(k; M) = P+(M) £fe , Ci u) (k; M) = P~(M)e k: (38) 

and the quasiparticle linewidth is in the form 

T«(M) = P-(M)W, 4 u) (M) = P+{M)W. (39) 

Here, P are the function of spin polarization 

^W-^- (40) 

In Fig. [1] we schematically show the spectral function for the up-spin electron A^(k,ou) = 
—lmG^(k,uj)/n calculated from eq. (|3~6|) . For the down-spin electrons, A^(k,uj) is obtained by 

replacing upper and lower subbands, e.g. zf = zi u \ etc. In the paramagnetic phase (Fig. 0(a)), 
quasiparticle dispersion are split into lower and upper subbands at u> ~ ±Jh with quasiparticle 
weight Z-r — Z<r = 1/2. As magnetization is increased below T c , the lower subband gains 

quasiparticle weight = (1 + M)/2 and that for upper subband decreases as = (1 — M)/2. 
At the ground state with perfect spin polarization M — 1, the electronic Hamiltonian describes a 
free electron system under Zeeman splitting field. Hence there exists only lower subband for the 

up spin electron. This limit is described £LS Z^ ^ — 1 and = 0. 

In the paramagnetic phase M = we see Xm S = — W , which means that the quasi-particle 
excitation is very incoherent; the lifetime of a quasi-particle is comparable with the time scale 
that an electron transfers from site to site. This result justifies us to take the D = oo limit which 
is essentially a single-site treatment. 

In the following section ( §3.1| ), we will discuss the behavior of the spectral function for general 
cases. 

Infinite Jh limit —Curie temperature—. In the semicircular DOS case, we consider the 
ferromagnetic state under doping at Jh 3> W. We set 

CfT 1 (iw n ) = {}UOn + Rn)I + Qn^z (41) 

and fi = — Jh + In order to keep the carrier concentration finite, we take the limit Jh — > oo 
with keeping Sfj, = 0(W). In this limit, we have 

Rn 8 \z n + R n + m z Q n ) ' ® n 8 \z n + R n + m z Q n ) ' 



where z n = iu n + Sfi. The Boltzmann weight is calculated from eq. ( |T2"D as 

(m) cx exp 



E 1 ^ 1 



(43) 



This equation tells us that the model is not simply mapped to the Heisenberg model. The 
Boltzmann weight of the Heisenberg model in infinite dimension is expressed as Pneis(^) ~ 
exp(— flh e gm z ), which contradicts with that of the DE model, i.e. P(m) ^ Pn C i S - This reflects 
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0.06 




Figure 2: Curie temperature for the semicircular DOS at Jh = oo (solid curve). Dotted 
curve is the form T c oc x(l — x). 



the fact that the itinerant ferromagnet has smaller T c due to the spin fluctuation affecting the 
itinerant electrons, compared to the insulating ferromagnet with the same spin stiffness. 
Let us calculate T c . At T ~ T c , we have M = (m) <C 1 and 



By solving eqs. 



Qn 

and 



W 2 



8 z n + R n 
W 2 



O(M), 



m 2 



z n -\- R r 



\Zn Rn 



;Q n )+0(M 2 ). 



(44) 
(45) 



we have 



Rn = V^n 2 - W 2 /A - Z n , 



Qn 



Rr, 



M 1 -8(m z 2 )Rl /W 2 ' 
Here, (m z 2 ) = J dQm z 2 = 1/3 at S = oo. The Boltzmann weight is given by 

Q, 



(46) 



P{m) oc exp 



E 



z n -\- R r 



-m. 



where 



J, 



ci-r 



■■ exp(-(3J eS Mm z 
8R 2 JW 2 



8R 2 J(3W 2 Y 



(47) 



(4J 



Then, the partition function is identical to that of the Heisenberg model with exchange coupling 
J e ff, and T c is obtained from 



f3^3W 2 - 8R 2 n 



(49) 



P=l/T c 



Solving eqs. (|46| ) and ( |49| ) self-consistently, we obtain T c as a function of /i, while the carrier 
number is calculated directly from Green's function G(iu} n ). In Fig. |2| we plot T c as a function of 
doping x. We see that T c has a maximum at x = 0.5 with T c ~ 0.05VT. Then, using (3W C 1 wc 
may approximate the summation in eq. (Ef9|) by integration, 



~ , 8R 2 (ix + Sa) 
dx- 



2tt J-oc 3W 2 - 8R 2 (ix + 6fi) 



(50) 



where R(z) 



z 2 -W 2 /4-z. 
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The result is particle-hole symmetric, namely T c (x) — T c (l — x). At x — > and x — > 1, T c 
diminishes because the ferromagnetism of this system is due to the kinetic energy of the conduction 
electron. In Fig. |2|, we also depict a curve T c oc x(l — x) proposed in ref. [34] for comparison. This 
simple form roughly reproduces the doping dependence of T c . In the following section, we discuss 
more precisely about the Curie temperature at finite J-r/W. 



2.4 Monte Carlo method for finite size clusters 



On a finite-size clusters, it is possible to investigate the double-exchange system by numerical 
methods. The result is unbiased and exact within the numerical errors. 

The partition function of the present model with localized spins treated as classical rotators is 
defined by 

Z = Tr s Tr F exp (-(3[H({™i}) ~ ^}) , (51) 

where Trg and Tip represent traces over spin and fermion degrees of freedom, respectively. In the 
finite size system, Z is obtained by taking the trace over fermion degrees of freedom first and spin 
degrees of freedom afterwards. Fermion trace is directly calculated from the diagonalization of 
2N x 2N Hamiltonian matrix, where N is the number of sites. Trace over spin degrees of freedom 
is replaced by the Monte Carlo summation over spin configurations {fhi}. 

For a fixed configuration of classical spins {fhi}, the Hamiltonian is numerically diagonalized 
and we obtain 

2N 

Tr F exp(-/?[H({m i }) - fiN}) = Ul 1 + exp(-/3(^({m,,}) - //))] . (52) 

u=l 

E v (v — 1, . . . , 2N) are eigenvalues of the Hamiltonian matrix for a given configuration {fhi}. We 
have the effective action for the classical spin system 

S«({rht}) = - £ log (l + e -«^) , (53) 

V 

which gives 

Z = Tr s exp(-S cff ). (54) 
Monte Carlo update of spin configurations is performed using the Boltzmann weight of the state 

P({rhi}) oc exp(-S' eff ({m i })). (55) 

In a Monte Carlo unit step, orientations of each spins are updated using the Metropolis algorithm. 
Since the spins are classical, spin updates can be performed ergodically. 

Thermodynamic quantities are stochastically calculated. Quantities which are associated with 
localized spins are obtained directly from the thermal average of spin configurations. Electronic 
quantities are calculated from the eigenvalues and eigenf unctions of 7i{{fhi}). 

One of the advantages for taking the classical spin limit is that there exists no "negative 
sign problem", which is present in quantum spin models. In the classical spin limit, the spin 
degrees of freedom is completely decoupled from those of fermions, and the fermionic trace in 
eq. (|52"D is obtained by solving the noninteracting lattice fermion system with random static 
potential. Another advantage is that the real frequency dynamics of electronic properties are 
directly obtained, since eigenvalues of the Hamiltonian are calculated through the Monte Carlo 
procedures. There is no difficulties of analytical continuations form imaginary frequencies, as 
is present in some quantum Monte Carlo methods. A disadvantage of this method is that the 
auxiliary field is static i.e. non-local in imaginary time, so it is not possible to make^ local spin 
flip using the imaginary time Green's function as in the case for the Hubbard model.E3t§) 



-10- 




Figure 3: Monte Carlo results for DOS on 6 x 4 x 4 cubic lattice at Jh/W = 4, at f3 = 
W/T = 6 and 36. Error bars are within the symbol size. Curves in the figure show the 
DMF results, (left) Two peak structure at uj ~ ±Jh is seen for both DMF and Monte Carlo 
results, (right) Lower subband at u ~ — Jh- 



Let us now compare the dynamical mean-field theory and the cluster Monte Carlo method. 
Here we calculate the electron DOS 

= -^E JmG -(^ + ir ?)A (56) 

/co- 
on a finite size cluster system in Fig. |3|. We treat iV = 6 x 4 x 4 cubic lattice at Jn/W = 4 
and \i = — Jh ; where the dynamical mean-field approach gives T c = 0.028PF. In order to avoid 
delta-function singularities in finite size systems, we use an adiabatic factor rj = 10~ 2 to smooth 
the spectra. In the Monte Carlo calculation, we calculate at /3 = W/T = 6 and 36 , and compare 
with the dynamical mean- field result in the paramagnetic phase T > T c . The result shows that 
the DMF calculation is very accurate. 



III. RESULTS 

In this section, we show the results for the classical spin limit S = oo of the DE system using 
the DMF approach as well as the Monte Carlo calculation. Hereafter, the electronic bandwidth is 
taken to be W = 1 as a unit of energy. For the carrier electron number, we express by x = 1 — (n) . 
We define M = (m z ), where < M < 1. Also, we describe the total moment, or sum of moments 

of localized spin and electron spin, as M tot = (fw 2 + \<J^j- In a normalized form we describe 

M* = M tot /M sat , where M sat is the saturation value of M tot at the ground state. We make 
distinctions between M and M* , since analytical calculation is better understood by M while the 
comparison with experiments should be done by M*. Nevertheless, in the strong coupling region 
Jh 3> W we have M ~ M* so effectively there exists no major differences. 

3.1 Electronic structures 

Spectral function and the density of states. In Fig. ^ we show the spectral function 
for the up-spin electron A^(k,u) on a cubic lattice, where k/ir = (C? C? • From the particle-hole 
symmetry and the spin symmetry, down-spin part Ai(k,u) is reproduced by the relation 



A l (k,u)=A< [ (Q-k,-u), 



(57) 



A(k, co) 




co/W co/W 

Figure 4: Spectral function A-t(k,u) on a cubic lattice at Jk/W = 2 and x=0.3, (left) in 
the paramagnetic phase T = 1.05T C , and (right) in the ferromagnetic phase T = 0.5T C . Here, 

fc/T = (C,C,0- 



where Q = (tt, it, it). 

In Fig. £| we show the quasiparticle excitation structure and its temperature dependence. There 
exists two-peak structure at around u ~ ±J H . Above T c , peaks at upper and lower bands are 
symmetric and equally weighted. Below T c , the structure remains split but becomes asymmetric. 
For the up-spin electron, the integrated weight is transferred from upper band to lower band. We 
also see the change of the quasiparticle linewidth T. The lower band peak becomes sharper, which 
means the reduction of F or the enlonged quasiparticle lifetime. On the other hand, T for the 
upper band peak increases. 

To investigate the change of the spectral weight in further detail, we calculate the DOS A a (uj) 
by ^-integrating the spectral weight. In Fig. ||] we show the DOS as a function of temperature 
in the paramagnetic and ferromagnetic phases. Two subband structure at u ~ ±Jh reflects the 
quasiparticle structure as described above. At the ground state, lower subband is composed of 
up-spin only, and down-spin band exists only at the high-energy region. The bandwidth of the 
DOS becomes narrower as temperature becomes higher. The band center is fixed at ui ~ ±Jh 
which is the energy level of the atomic limit t = 0. 

These results are in agreement with the exact result in a limiting case (Lorentzian DOS with 
Jh — > oo) discussed in § |2.3| . Let us focus on the lower subband. The DOS is nearly proportional 
to the spin polarization of the localized spin, namely 

A^(uj) ~zf{M) =(l + M)/2, 

M u ) ~^°(M) =(l-M)/2, (58) 

where z is the quasiparticle weight discussed in § |2.3| . We see that z is determined by the magne- 
tization through (1 ±M)/2, which is the probability that the local spin is parallel to the itinerant 
electron with up (down) spin. This is easily understood from the nature of the double-exchange 
interaction which projects out the antiparallel component of the spins. This spin- dependent pro- 
jection may be viewed as the evolution of the majority-minority band structure. The up (down) 
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Figure 5: Temperature dependence of DOS for Ju/W = 2 and x = 0.3, where T c = 0.019VF. 
Peak structures are observed at oj ~ ±Jh. 



spin band becomes a majority (minority) band below T c or under magnetic field. Shift of the 
spectral weight from the minority band to the majority band occurs as temperature is lowered. 
At the ground state (M = 1), minority band completely loses its weight. 

Change of the electron bandwidth is understood qualitatively through Anderson-Hasegawa's 
picture.^ Electron hopping amplitude is proportional to cos(#/2) where 9 is the relative angle of 
the localized spins. At high temperature, 9 deviates from zero due to spin fluctuation, and the 
amplitude of electron hopping matrix element and hence the bandwidth decreases. This is also 
shown by the virtual crystal approximation.^ 

In order to account for the width of A(k,u), or the quasiparticle lifetime, we have to go 
beyond a mean-field picture like in Anderson- Hasegawa approach. The origin of the linewidth T 
is the thermal fluctuation of the spins. In the strong coupling region Jh ^> W, spin scattering 
phenomena in the paramagnetic phase is so large that quasi-particles lose their coherence due 
to the inelastic scattering by thermally fluctuating spins. In the ferromagnetic phase, the spin 
fluctuation decreases as temperature is decreased. For the majority band, this decreases T, and in 
the limit T —>■ the majority band becomes a free electron band. However, for the minority band, 
the spin projection causes further loss of the coherence, which leads to the increase of T as well 
as the decrease of the quasiparticle weight z. Asymptotically in the limit T — > 0, T approaches to 
a constant ~ W, and at the same time z ^ 0. Then, at T = 0, the minority band with finite T is 
projected out. 

Half metal. Metal with a DOS structureshown in Fig. [| where only one of the spin species 
have the Fermi surface is called a half metal£§EP Namely, because of the 'Zeeman splitting' due to 
Strong Hund's coupling Jh 3> W, ferromagnetic ground state of the DE model shows a perfect spin 
polarization and thus is a half-metal. Experimentally, spin-resolved photoemission investigation 
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00) shows that the conduction band of the doped manganites is a half- metal. Artificial trilayer 
junction of manganitesB' also shows a large tunneling magnetoresistance phenomena, and the 
spin polarization is estimated to be more than 80%. Such a DOS structure creates a phenomena 



called tunneling magnetoresistance (TMR), which will be discussed in §5.1 



Shift of the chemical potential. A direct consequence of the change of the bandwidth 
controlled by the magnetization will be observed in the shift of the chemical potential. In this case, 
the change of the DOS structure is in a way such that the band center is pinned by the Hund's 
coupling energy ± Jh and the band edge shifts away from the center as magnetization is increased. 
Then, for a hole doped case the position of the chemical potential increases by increasing the 
magnetization. 
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Figure 6: Temperature dependence of fx at Jh/W = 4 and x = 0.20 under various magnetic 
field. Inset: Afi/W as a function of M 2 . Lines show the result at T < T c for H = 0. Squares 
and diamonds are data at T = 1.01T C and 1.2T C by applying H, respectively. 



In Fig. |^ we show the temperature dependence of the chemical potential at x = 0.2 under 
various magnetic field-E^P At H = 0, chemical potential fx is nearly temperature independent 
above T c . Below T c , fx shifts as a function of temperature. We also calculate fx and M, (i) at 
H = by changing temperature in the region T < T c , and (ii) at fixed temperature above T c by 
changing H. In the inset of Fig. ^| we plot fx as a function of magnetic moment M 2 for both cases. 
As a result, we see the scaling relation 

Afi/W oc M 2 , (59) 

where A fx = [x(T, H) — fi{T = T c , H = 0). We see that Afx can be as large as 0.1IT. 

Thus, for a fixed band filling, the total change of the DOS width in the entire energy range 
causes the shift of fx. The change in such a large energy scale controlled by magnetization produce 
the characteristic feature of the shift of fx in DE systems; namely, that the shift of fx is as large as 
a few tenth of W and that the scaling relation (|59|) is satisfied up to such a large energy scale. 

Such a large shift of fx might possibly be applied to electronic devices which controls the MOS 
gate voltage by the magnetic field. 

3.2 Magnetic structure and transport properties 

Magnetic transition temperature. In the limit Jh — > oo, Curie temperature T c of the 
DE model is determined by the electron kinetic energy. Indeed, DMF calculation shows that T c 
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i.e. 



T c oc W for Jh = oo limit (see §2.3). Here we show the case for 



is scaled by electron hopping, 
finite J H /W. 

In Fig. |7| we show the Curie temperature T c as a function of doping x for various values of 
Jn/W. At finite Jn/W, T c is reduced from Jh = oo values. We also see that T c systematically 
increase as x is increased and have maximum at around x ~ 0.5, which is due to the increase of 
the kinetic energy. 
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Figure 7: Curie temperature T c as a function of Jn/W and x. 



At half-filling x = 0, there exists antiferromagnetic order. In Fig. || we show the phase diagram 
at x = for the D = oo Bethe lattice. In the weak coupling region Jh <C W, the Neel temperature 
Tn is equivalent to the results from the SDW mean-field type equation with Jn/W dependence 
in an essential singular function. At Jh 3> W, Tn is determined from the Heisenberg model with 
the exchange coupling Jaf ~ t 2 / Jh oc W 2 / Jh- 
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Figure 8: Neel temperature at x = 0. 



Thus we see the antiferromagnetic order at x = and ferromagnetic ground state at sufficiently 
doped case. In the underdoped region x <C 1, de Gennes discussed the existence of the canted 
statetP. However, in §|3]4| we show that the canted state is unstable against phase separation, i.e. 
mixed phase of x = antiferromagnetic region and x > ferromagnetic region.Ej) 

Resistivity as a function of magnetization. Resistivity p(T) as well as total magneti- 
zation M tot as a function of temperature is given in Fig. ^[ Here, po is a constant of resistivity 
which corresponds to the Mott's limit value (inverse of the Mott's minimum conductivity) in three 
dimension. M sat is the saturated magnetization at T = 0. 



- 15 - 




Above T c , the value of the resistivity is in the order of the Mott limit p(T) ~ po, with small 
T- dependence- Below T c , resistivity drops quickly as magnetization increases. From the DMF 
calculation,&y) it has been made clear that the resistivity behaves as 

p(M)/p(M = 0) = 1 - CM 2 , (60) 

where C is a temperature/field independent constant. Namely, temperature and magnetic field 
dependences come from the magnetization M = M(H,T). As a function of magnetization M, 
all the T and H dependent values p(T, H) converge on a universal curve in (|60|) . In other words, 
the origin of the resistivity is due to spin fluctuation, oiLmore precisely spin disorder scattering, 
discussed by Kasuyac3) and later by Fisher and Langer.fcf 

In the Born approximation (weak coupling limit), we see C — A phenomenological treat- 
ment by Kubo and Ohata which estimates the resistivity from the spin fluctuation, 

p oc (5S) 2 oc 1 — M 2 , (61) 

also gives C — 1. However, in the DMF at Jh 3> W, we have C > 1 which indicates the strong 
coupling behavior. In the next section, the relation with experimental MR is discussed in more 
detail. 

Above T c , small T-dependence in p(T) is observed within the DMF treatment. This result 
might be an artifact of the approximation since the local spin fluctuation is saturated above T c . 



3.3 Charge and spin dynamics 



Optical conductivity. Temperature dependence of the optical conductivity is shown in 
Fig. [10]. In the paramagnetic phase, the spectrum splits into two peaks due to the 2-subband 
structure of the DOS. Namely, intraband particle-hole channel creates a Drude-like peak at uj ~ 0, 



while interband channel creates a peak at around oj ~ 2 Jh- In the inset of Fig. [10], we show the 
weight of the interband process as a function of 1 — M* 2 . The integrated weight of the interband 
optical process at u ~ 2Jh defined by 

S= du <t(u), (62) 

where cutoff frequency is taken as u> c = Jr. We see a scaling relation 

S oc 1 - M* 2 . (63) 
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Figure 10: Temperature dependence of the optical conductivity for Jn/W = 2 and x = 0.3. 
Inset: Integrated weight S at high energy part uj ~ 2Jh- 



This is explained by the temperature dependence of DOS. Below Curie temperature, the DOS 
changes as in Fig. |5|. Interband optical process at w ~ 2Jh is constructed from a process of 
making a pair of lowerband hole and upperband electron. Then, the transfer of spectral weight 
by magnetization creates the change in the optical spectra as follows. For the up spin electrons, 
the spectral weight of lower and upper subbands are proportional to (1 + M*)/2 and (1 — M*)/2, 
respectively. Since this optical process conserves quasiparticle spin, the total weight of the optical 
spectrum is proportional to the product of the initial-state weight at lower subband and the final- 
state weight at the upper subband, and hence scaled as 1 — M* 2 . Contribution from the down-spin 
band is also the same, and we have S oc 1 — M* 2 . 

Stoner excitation. Stoner susceptibility is calculated by 

z) = — G ( k + Q, i^n + z)G(k, iuj n ). (64) 

Here, correlation effects are taken into account through the self-energy correction in G. In Fig. [ll] 
we show g-dependence of lmx(q,^>) for various temperatures, at Jr/W = 2 and x = 0.3. We see 
two-peak structure at uj ~ and uj ~ 2 J H which is explained from the Jn-split DOS structure. The 
Stoner absorption is produced from a particle-hole pair excitations with spin flip, which produces 
a peak at low energy from intra-band processes and another peak at uj ~ 2Jh from interband 
processes. 

We see a weak q dependence in the low frequecy part, especially at T T c . This part of the 
Stoner process is dominated by the combination of the majority- minority quasiparticles. Since the 
quasiparticles in the minority band is incoherent, Imx{q,u) is weakly q dependent. On the other 
hand, high energy part of Imx{q, cj) at T <C T c have larger q dependence. This part is dominated 
by the majority- majority quasiparticle channel. Thus q dependence of Xmx{q, reflects the band 
structure of the quasiparticles. 

Let us see the u dependence at the low frequency region. In Fig. [12] we show Xzb(^) = x{Qi u ) 
where Q = (tt, 7r, tt). We see that at small uj we have tu-linear relation, i.e. Tmx cf-UJ at uj W. 
Coefficients for a;-linear part decrease by decreasing the temperature, and we findfi^ 

Imx(Q,to) oc (1 - M* 2 )uj (65) 

for small values of uj. The relation ( |S5D is observed at all values of q with weak q dependence. 

Low energy Stoner absorption is constructed from a minority particle and majority hole chan- 
nel. Since the majority and the minority bands have the spectral weight proportional to (l+M*)/2 
and (1 — M*)/2, respectively, the low energy part of the Stoner absorption is proportional to their 



-17- 



ffl/W T=1.05T C T=0.25T C 




rxMRrrxMRr 



Figure 11: Contour plot of the Stoner absorption x(?j w ) on a cubic lattice at Jn/W = 2 
and x = 0.3, in the paramagnetic phase T = 1.05T C (left) and in the ferromagnetic phase 
T = 0.25T C (right). At T = 0, the low energy part at u> ~ disappears. 
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Figure 12: Stoner absorption x(Q,u;) at the zone boundary Q = (7r,7r,7r). At Jn/W = 2 
and x = 0.3, transition is at T c /W = 0.019. 



product, 1 — M* 2 . The cj-linear behavior comes from the Fermi distribution function. The in- 
coherence of the minority band gives the weak q dependence. Thus we have the scaling relation 
ImX oc (1 — M* 2 )uj. The weak q dependence is in large contrast with the conventional weak 
ferromagnet where minority band is also coherent, which gives strong q dependence through its 
band structure. 

3.4 Phase separation^ 

Magnetic phase diagram of the weakly doped DE model has been studied by de Gennes.i- 1 
Assuming the homogeneity of the doped carriers, he concluded that the spin canted phase is the 
most energetically favorable state. However, we have recently shown that there exists an instability 
toward phase separation,!!^ and the assumption of uniformly doped charges by de Gennes is not 
valid. 

One of the ways to discuss the phase separation is to make a grand canonical calculation of the 
particle number x as a function of the chemical potential fi. If there exists a jump of x(/a) at the 
critical value fi = /i c in the thermodynamic limit, it implies that two phases with different doping 
x coexist at /i = /i c . We have showrM the jump of x(/i) in the DE model for sufficiently large 

*Works shown here concerning the issue of phase separation has been made in collaboration with S. Yunoki, A. 
Moreo and E. Dagotto. 
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Jn/W. The jump occurs from x = to a state with finite x. The calculation has been performed 
within the DMF approach (D = oo), and cross checked by the Monte Carlo method for the D = 1 
and 2 clusters in the extrapolated limit of T — > 0. 

Phase boundary is simply determined from the jump of x(fi). 
phase diagram. 



In Fig. |3 

At the low temperature region, we see mixed phases of x 



we show the x-T 



AF state and doped 

(x > 0) state with either paramagnetic or ferromagnetic state, depending on temperature. 




Figure 13: Phase diagram at J-r/W = 4. P (F) represents paramagnetic (ferromagnetic) 
region. At x = we have antiferromagnetic (AF) phase. Regions labeled by AF-P (AF-F) 
are phase separated region with mixed phases of AF and P (AF and F) phases. 



The mechanism of the phase separation may be understood as follows. In Fig. [14] we show 
the change of the density of states for different doping concentration. We see that the bandwidth 
substantially differs for different magnetic states, while the band center remains at u ~ ±Jh- As 
discussed in §3.1| , this typical DOS structure is a consequence of the DE half-metallic system. Let 
us consider the zero temperature limit. At the x = AF state, we have \l = 0. In order to hole 
dope this AF state, we need to decrease /j to a AF gap- value, /i = — Aaf- However, as we see 



in Fig. [14], it is also possible to make a F state at fi = —Aaf with rich hole density x. AF state 
gains the exchange energy ~ t 2 /Jn, while the F state gains the kinetic energy ~ tx. Thus, in the 
limit Jh ^ t, the doped F state becomes energetically favorable. In such a case, there exists the 
chemical potential /i c such that —Aaf < < 0, where the energies of AF state with x = and F 
state with x > are equal. Level crossing from AF to F states occurs at this critical point /i = fi c . 
Doped F state is realized before doping the AF state. Thus we have a jump in the hole density 
from x = to finite x. 

The discussion based on the DOS structure is quite generic. In strongly correlated systems 
with bistable phases, macroscopic change of the order parameters creates a large change in the 
electronic DOS. This gives a macroscopic jump in (n) from one phase to another when fi is fixed. 
Then, phase separation is associated with the density- driven phase transition. 

Let us discuss the issue from a different viewpoint. In Fig. [15] we show the equal time spin 
correlation 

%) = ^E(^-^)e i9(i - j) (66) 

ij 

on a one dimensional system at L — 40, Jb_/W = 4 and (3 = 150iy _1 . We clearly see the 
crossover from the antiferromagnetic state at n = 1 with a peak of S(q) at q — 7T, and the 
ferromagnetic state at n ~ 0.7 with the peak at q — 0. Around n ~ 0.9, where large change of 
n(/j,) is observed, we see two peak structure at q = ir and q = 0. This result is consistent with 
mixed phase of ferromagnetic and antiferromagnetic states. On the other hand, magnetic states 
with incommensurate momentum, typically at q = 2kp, have been observed for the weak coupling 
region Jh ^ W. Thus the tendency toward phase separation is prominent at the strong coupling 
region. 



- 19 - 





.". 




1 li 1 




; ; 


--AF 




"c 








3 


; ! 


— F 




(arb 


/ i 




DOS 




. A : 


ra 


; A 




To! 







-4-3-2-101234 

co/ W 



Figure 14: Density of states for ferromagnetic (F), paramagnetic (P) and antiferromagnetic 
(AF) states at x = 0.2, 0.1 and 0, respectively, for J H /W = 2 and T/W = 0.005. 




Recently, the issue of phase separation is also investigated by_yarious methods.c3'li§E!P The 
model in the weak coupling region Jh <C W has also been studied.B* However, the mechanism of 
the phase separation might be different from the strong coupling region where the half-metallic 
DOS plays^nimportant role. In manganites, several experiments claim the existence of phase 
separation.E3&P We will later discuss this issue in comparison with experiments in detail. 



IV. 1/S CORRECTIONS 

Spin wave expansion. We introduce a linear spin wave theory. The spin wave operators 
are introduced from 

Si + ~ V2Sa i} Sr ~ V2Sa\, S* = S - a}^. (67) 

Hereafter we restrict ourselves to the lowest order terms of the 1/S" expansion at T = 0, where the 
localized spins are perfectly polarized so that (aja.j) = 0. We consider the half-metallic ground 
state, i.e. fki = where fk a is the Fermi distribution function. 
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The spin wave self-energy in the lowest order of 1/S expansion is obtained diagramatically 
as§) 
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(61 



Here, G a (k, iu n ) = {iuj n — e k + cr Jh) 1 is the fermion Green's function, and ft = 1/T. 

The spin wave dispersion relation oj q is obtained self-consistent ly from oj q = H(q,uj q ). We have 
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(69) 



where is the fermi distribution function of the majority band. 

Let us consider the strong Hund coupling limit Jh 3> t. If we assume a simple cubic lattice 
with nearest-neighbor electron hopping, 



we have 



where E* 



2t (cos k x + cos k y + cos k z ) , 

3 — cos q x — cos q y — cos q z 
6 



u q= Q — oj q= o is the spin wave bandwidth, given by 
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(70) 
(71) 

(72) 



We see that in the strong coupling region the spin wave bandwidth is determined only by the 
electron transfer energy. 

In the isotropic case, the spin stiffness is defined via uo q = Dq 2 in the long wavelength limit 
q — ► 0. From eq. (|69|) we have 



D 



1 1 

2SN 



ld 2 e k 



1 



' de k 



2 «% 2 2J H V ^ 



(73) 



We emphasize that the expansion given here is not with respect to Ju/(tS) but to 1/S, due 
to energy denominator 2Jh between up- and down- spin electrons. This is understood from the 
fact that Jh plays an role of a projection and does not enter the energy scale by itself in the limit 
Jh ~~ > 00. Therefore, the calculation is valid even in the large Hund's coupling limit Jh ^> t as 
long as we restrict ourselves to T = 0. Indeed, the result for J H — > 00 obtained by the present 
approach is equal to those in the projection limit Jh = 00 shown by Kubo and OhataHD 



Discussion. In the DE model with Jh 3> W we observe a short-range spin interaction, in 
contradiction with the case of Jh -C W (s-d model) where the well-known RKKY interaction is 
long ranged with a power-law decay. The difference comes from the electronic structure. 

In the case Jh <C W where electronic DOS is the ordinary one (not half-metal), the gapless 
quasi-particle excitation of the electronic part creates a particle-hole spin excitation channel with 
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2/cf singularities. Interactions among localized spins are mediated by this gapless channel and 
thus long ranged. 

On the other hand, for the DE model with Jh 3> W the half-metallic DOS creates a gap due to 
Jh splitting, and unlike the former case the particle-hole spin channel is massive and short ranged. 
The qualitative difference comes from the half-metallic structure. In the real-space picture, we 
consider a perfectly spin polarized state at T = and twist a spin at site iq. Spin polarization of 
itinerant electrons are along the total polarization axis except for the ig-th site where it orients 
toward the local spin direction S io . In the strong coupling limit Jh W, the electron at site i 
is localized because it has different spin orientation from the spins in neighboring sites. Since the 
effective interaction between localized spins are mediated by the motion of electrons, the effective 
spin-spin interaction is short ranged. As Ju/W increases, electrons become more localized so the 
range of effective interaction becomes shorter. In the extreme limit Jh — > oo, the interaction is 
nonzero only for the nearest neighbors which gives a cosine-band dispersion. 

Let us discuss the higher orders of 1/S expansion terms. In an ordinary (insulating) Heisenberg 
ferromagnet, the first 1/S term is the relevant term with respect to the one-magnon dispersion uj q . 
Higher orders of 1/S* expansion only give magnon-magnon interaction terms, and thus irrelevant 
within the one magnon Hilbert space. However, for the DE model, higher order terms also give 
one-magnon kinetic terms. We need to take into account the asymptotic 1/5* expansions even 
for the one-magnon dispersion relations. Such higher order terms may be considered as vertex 
corrections to the self-energy term H(q, u). 

Numerically, Kaplan et alB studied the S = 1/2 case at Jh ^ oo. They observed cosine-band 
type behaviors in the well-doped cases, and the deviation from them in the limit n — > and n — *■ 1. 
The result might be understood from the Migdal's discussion. Let us consider the electron kinetic 
energy E^ n and the energy scale of magnons (u q ) of the DE model. The deviation from a cosine- 
band comes from the vertex correction which is relevant if ^ which occurs in the lightly 
hole/electron doped region n — * 1 or n — > 0. On the other hand, in the well-doped region we have 
-Ekm r<S (uj q ) and the vertex corrections are small. In the region where doped manganites show 
ferromagnetism, vertex corrections do not seem to be important. This explains the consistency 
between experiments and 1/S" results 



V. COMPARISON WITH EXPERIMENTS 

Abbreviations: Hereafter we use these abbreviations: 
La1_a.Sra.MnO3 (LSMO), Lai-aPb^MnOg (LPMO), La^Ca^MnOg (LCMO), 
Pr^S^MnOs (PSMO), Nd^Sr^MnOa (NSMO), Pn_a.Ca_M_.O3 (PCMO). 

5.1 Experimental —varieties of properties in "CMR manganites"— 

Let us briefly mention the varieties of phenomena in CMR manganites. For a review of recent 
experiments, readers are referred to ref. 0. It is emphasized that systematic studies of A-site 
substitution is quite important to understand the complex behaviors in manganites. Extrinsic 
effects due to grain/domain boundaries are also discussed. 

A-site substitution. Recent improvements in precise control of the A-site cations substitu- 
tions in AMnOs revealed a complex phase diagram as a function of substitution, temperature and 
magnetic field. They exhibit various phases with magnetic, charge, orbital and lattice orderings. 
For example, phase diagram for doping (x) vs. temperature (T) is well known for Lai_a;Caa;Mn0 3 
(LCMO),!) as well as L_1__Sr._M_.O3 (LSMO)|§ Ndi__Sr_Mn0 3 (NSMOLand Pri_ a Ca_M_0 3 
(PCMO).S Effects of A-site substitution is also studied for a fixed doping.BSi^ 

Major effects of A-site substitution are the bandwidth control and the carrier number control. 
It is well understood that the ratio of rare-earth (3+) ions and alkaline-earth (2+) ions determines 
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(a) (b) 

Figure 16: (a) Schematic phase diagram at x ~ 1/3 by ionic radius of A-site {r\). Abbre- 
viations are: paramagnetic metal (PM), paramagnetic insulator (PI), ferromagnetic metal 
(FM), charge ordered insulator (COI), as well as Curie temperature (T c ) and charge or- 
dering temperature (Too)- (b) Qualitative behaviors in p{T) for (La,Sr)Mn03 (LSMO), 
(La,Ca)Mn0 3 (LCMO) and (Pr,Ca)Mn0 3 (PCMO) at x ~ 1/3. Uparrows (f) in the figure 
show T c , while downarrows (!) indicates Tqq. 



the nominal values of the carrier number x. At the same time, change of the average radius of 
the A-site ions (r\) by chemical substitutions gives the "bandwidth control" through the chemical 
pressure. Such kind of chemical control creates a largechange in the nature of the compounds. 
In general, compounds with larger (r^) have higher T c .E3'l3'EZP It is considered to be due to wider 
effective bandwidth for e g electrons in larger (r\) compounds, since it gives less Mn-0 octahedra 
tilting. 

However, we should note that it is still controversial whether the phase diagram is controlled 
mostly by the bandwidth alone. For example, ionic size variation cr(rA) also plays some role to 
change T^H* as well as the fact that decrease of T c for small (ta) is substantially larger than the 
estimated value from the change of bond angles. 

Let us concentrate on the region x ~ 1/3 where it is far from antiferromagnetic insulating 
phase at x ~ and the region with charge and orbital ordering at x ~ 0.5. The compounds are 
roughly classified as follows: 

• High T c compounds: e.g. LSMO. 

A canonical example for the high T c compounds is (La,Sr)Mn03 (LSMO) with T c ~ 380K. 
Resistivity shows a small value at lowest temperature (p ~ 10 2 /if2cm). At T ~ T c , p{T) 
takes much larger value but still in the order of Mott's limit p = 2 ~ 4mf2cm. Above T c , 
p(T) shows a metallic behavior, i.e. dp(T)/dT > 0. Namely, this compound is a good metal 
below T c and become an incoherent metal above T c with the absolute value for p(T) being 



near Mott's limit. 
Low T c compounds: e.g. LCMO, PSMO. 

LCMO is the most well-investigated compound. Ca substitution creates smaller (ta) and 
larger cr(rA)- It has lower 280iT compared to LSMO, and shows metal to insulator 

transition at around T c . 

Compounds with charge ordering instability: e.g. PCMO. 
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As (ta) is further decreased, compounds show a charge ordering at T ~ 2OOK.i30il'li ) In 
the zero-field cooling process, ferromagnetic metal phase does not appear. 

At x ~ 1/2, the phase diagram becomes more complex. At the lowest temperature, the 
tendency of the competition between ferromagnetism and charge ordering!!' © driven by (ta) 
control remains the same. It has been recently discovered that in the intermediate region, d^jiew 
phase of A-type antiferromagnetic metal region exists in the narrow vicinity of x = 1/2.E3EP 
Temperature dependence is also complex. Behaviors of the resistivity above T c are roughly the 
same with the case of x ~ 1/3. For example, metallic behavior above T c is observed for NSMO 
at x = 1/2 as well as compounds with larger (ta)- Substitution of Nd by Sm reduces (ta) and 
makes the paramagnetic phase insulating.^ 



Extrinsic effects in polycrystal samples. In samples with multiple grain structures, it 
has been shown that there £xist so-called tunneling magnetoresistance (TMR) phenomena through 
the spin valve mechanism.^'^tBEP Magnetoresistance of the material with artificially contrplled 
grain/interface boundaries are also studied to realize a low-field MR device through TMR.c30E3) 

In manganites, the half-metallic behavior due to DE interaction is considered to cause a large 
amplitude in MR. Spin polarizations in grains are schematically depicted in Fig. [17|. In the zero 
field case, each grain have different spin orientations. Namely, the spins of the half-metallic 
electrons are different from one grain to another. Due to the nearly perfect polarization nature of 
itinerant electrons, inter-grain hopping amplitude is suppressed by such random spin polarizations. 
Under the magnetic field, spin polarization of grains become parallel to the external field. Inter- 
grain electron hopping becomes larger in this case. For multi-grain system with half-metallic 
states such spin valve phenomena becomes prominent and gives MR effect in a low field range. 

Resistivity in polycrystal samples seems to be dominated by such extrinsic effects. One should 
be careful about discussing the experimental data from a microscopic point of view. Such TMR 
behavior in these perovskite manganiles_should be discussed in relation with other half-metallic 
materials such as Cr02 or Tl 2 Mn 2 O7.S0 





H = 



H>0 



Figure 17: Grains of the polycrystal samples, (a) at zero magnetic field H = 0, (b) under 
magnetic field H > 0. Arrows indicate the orientation of the magnetization for each grains. 
For the relation with low field TMR, see the text. 



5.2 Comparison with high Curie temperature compounds 

Here, we will show the comparison of theoretical results with experimental data of high T c 
compounds such as LSMO and LPMO. We discuss that the DE Hamiltonian alone explains most 
of the thermodynamics of these manganites, including Curie temperature and resistivity. 
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Curie temperature. In Fig. [18] we plot the x dependence of the Curie temperature for 
LSMO, together with the fitting curves obtained by DMF.tP From the fitting, we see that with 
parameters W ~ leV and Jn/W ~ 4 the value of T c as well as its x-dependence is reproduced. 
The bandwidth of W ~ leV is a typical value fiir^Sd-transition metal oxides, and is consistent with 
the band calculation estimate for juanganitesJZj'LJ'lZ§' It is also consistent with the value obtained 
from the spin wave dispersion fit,£3) shown later in this section. 
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Figure 18: Curie temperature T c of the DE model in comparison with those for 
Lai-xSr^MnOa. Experimental data are from ref. 77. 



Thus the value of T c in LSMO is reproduced by the DE model alone. This result is consistent 
with other methods, e.g. high-temperature expansion by Roder et al& as well as by the Monte 
Carlo method, yp 

Millis et alu discussed that the magnitude of T c as well as its x dependence for LSMO cannot 
be accounted for by the DE model alone, and discussed the importance of the dynamic Jahn- Teller 
effect. However, this part of their discussion is due to an inappropriate estimate of T c which is 
based on a calculation of effective spin coupling at T = 0. In the itinerant systems, T c is reduced 
by spin fluctuation, in general. Although T c scales with W, it is quite small compared to W due 
to the prefacter, whose typical value is T c ^ 0.05W 7 (See Fig. 0). Thus bandwidth of leV creates 
T c of the order of room temperature. Decrease of T c as x decreases is explained by the reduction 
of the kinetic energy!^ 



Resistivity and magnetoresistance. Resistivity of the high-T c compounds in single crys- 
tals at sufficiently large doping x ~ 1/3 are different from those of low-T c compounds or polycrystal 
samples. For LSMO at x ~ 1/3,© residual resistivity p is m the order of a few 10pf2cm, and 
the temperature dependence p(T) shows a monotonously increase, i.e. dp/dT > even above T c . 
(For LPMO, see ref. [79|. Recent experiment by Cheong et a/.© shows that the resistivity of LSMO 
continuously increase up to 1000K, without saturation or metal-semiconductor transition.) The 
value of resistivity at T c is typically p(T c ) = 2 ~ 4mOcm, which is in the order of the Mott limit. 
In short, LSMO is a good metal at T T c , and an incoherent metal at T ^ T c . The DE model 
reproduces these data (See Fig. |^). Similar temperature dependences of resistivity are observed 
in a wide class of materials of half-metals such as Cr02 and Heusler alloys.0- ) 

It has been discussed that the DE model cannot explain the resistivity of LSMO in its absolute 
value as well as the temperature dependence.© However, it is now clear that if one compares data 
for a high quality single crystal of LSMO (not polycrystal, or other compounds with lower T c ), 
DE alone does account for the resistivity. 

In Fig. [II] we show magnetoresistance of LSMO at x — 0.175. Universal behavior of the 
magnetoresistance in the form 

- Ap/po = CM* 2 (74) 
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Figure 19: Magnetoresistance, plotted by magnetization vs. resistivity. Solid curves and 
filled symbols are the data for Lao.875Sro.i7sMn03 from ref. JF?]. Open symbols and grey curve 
show DMF results for Ju/W = 4 and x = 0.175, for semicurcular (S-DOS) and Lorentzian 
(L-DOS). 



is observed, where C ~ 4 is temperature/field independent constant. US Namely, the resistivity 
is directly related to the magnetization. In Fig. |19] we also show the DMF results. The curve 
reproduces the experimental data.E-tf* The relation Ap/p(0) = —CM 2 is obtained theoretically in 
eq. (pop . The value C > 1 shows that the system is in the strong coupling region. 
Note that in LCMO polycrystal samples, a different scaling relation 

p(M) = p(0) exp(-aM 2 /T) (75) 

has been reported.® From such an activation-type temperature dependence, polaronic origin of 
the resistivity has also been discussed. The difference between the scaling relation in LSMO ( j74|) 
and LCMO (|75|) indicates the qualitative difference for the origin of CMR. 

Thermoelectric power. Seebeck coefficient S for LSMO has been reported-H* They show 
non-universal behaviors, including the change of the sign. However, in the vicinity of T c , a scaling 
behavior in the form —AS/ S(0) ~ — Ap/ p(0) irrespective of doping is reported, where AS and Ap 
are the change of Seebeck coefficient and resistivity under magnetic field, and S(0) and p(0) are 
their zero-field values, respectively© In Fig. |20]we plot the DMF result for the Seebeck coefficient 



S on a Lorentzian DOS at Jn/W = 4. The data is plotted in the form —AS/S(0) vs. — Ap/p(0). 
We see — ASyS'(O) ~ — Ap/p(0) for different values of x. 

Spin excitation. From the neutron inelastic scattering experiment, spin wave dispersion 
relation of Lao.7Pbo.3Mn03 (LPMO at x = 0.3) is investigated.^ LSMO in the ferromagnetic 



metal region has^also been studied.0' 

Perring et alx3 found that the experimental data of spin wave dispersion relation for LPMO 
fits the cosine band form. They phenomenologically argued that a ferromagnetic Heisenberg 
model with nearest-neighbor spin exchange couplings 2J ef[ S ~ 8.8meV is a good candidate for the 
effective spin Hamiltonian, although the_ma±ejnal is a metal. 

Now we discuss the DE model results.l3ErEj In Fig. |21, we plot the theoretical results together 



with the neutron inelastic scattering experiment data. As we have discussed previously, cosine- 
band dispersion is obtained in the limit Ju —>■ oo. This gives the identical fit with the analysis by 
Perring et ai, and the fitting parameter determines electron hopping as t ~ 0.3eV. Here, finite 



gap at q = is artificially introduced as done in ref. |82], due to some experimental inaccuracies. 

At the zone boundary q ~ (rr, rr, rr) we see the deviation from the cosine band, or softening 
of the dispersion. The data is well accounted for by introducing some finite value for Jn/t. In 
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Figure 20: Seebeck coefficient under magnetic field, scaled by resistivity. The line is a guide 
to eyes. 
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Figure 21: Spin wave dispersion at x = 0.3. The experimental data for Lao.7Pbo.3Mn03 is 
from ref. 



Fig. ^1] we show the result for Jn/t = 12. From fitting we obtain t ~ 0.26eV and Jh ~ 3.1eV. 
Although it is not easy to estimate the systematic errors of the fitting, it is plausible to say that 
the DE model with t = 0.2 ~ 0.3eV and Jh ^ 3eV explains the spin wave dispersion relation of 
LPMO. These values are the effective hopping energy of the double-exchange model with reduced 
degrees of freedom. —Nevertheless, these values of t are consistent with those estimated from the 
band calculation.L3'll3L3) 

Anomalous damping of zone boundary magnons are also reported.!^ Within the DE model, 



spin wave excitation at finite temperature interacts with the Stoner continuum as shown in Fig. |12 
From eq. (j65|) , we see that 



T(q,T) oc (1 -M 



*2~ 



(76) 



where T(q,T) is the linewidth (inverse lifetime) of the magnon. This explains that the magnons 
are damped at finite temperature, especially at the zone boundary. Indeed such a scaling relation 
([76} is observed in LSMOH 

Optical conductivity. Optical conductivity for Lao.eSrcuMnOs by Moritomo et al& is 
shown in Fig. ^(a). Here, temperature independent part which is discussed to be due to d-p 
charge-transfer type excitation is subtracted. There exists a peak at u> ~ 3eV, and its temperature 
dependence is in a way that it vanishes at T — > 0. 
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(a) Experiment La 6 Sr . 4 MnO3 r c =325 K 
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Figure 22: Optical conductivity, (a) Experimental data for Lao.6Sro.4Mn03 taken from 
ref. 87, where temperature independent part is subtracted as a background, (b) DMF result 
for W = leV, J H = 1.5eV and x = 0.4. 



By choosing W = leV, Jh = 1.5eV and x = 0.4, it is possible to reproduce the experimental 
data by the DE model. In Fig. ^(b), we show the DMF results for the semicircular DOS. The 
peak structure as well as its temperature dependence is in agreement. 

Let us discuss more details about the scaling relation.ll3'EZP Experimentally, the integrated 
spectral weight 

(77) 



<? = 

'-'exp — 



duja(uj) 



2.2eV 



shows a scaling relation 



S exp oc 1 - M* 2 . 



(78) 

This is explained by the transfer of DOS by magnetic fluctuation in the DE model, as shown in 
eq. (H). (See also Fig. in 

Thus, optical spectra show experimental evidence for the shift of DOS which is a typical 
phenomena in DE model, and give the estimate of the parameter Jh ~ 1.5eV. 



Summary: LSMO as a canonical DE system. As long as the high-T c compounds are 
concerned, e.g. LSMO at x — 0.2 ~ 0.4, the DE Hamiltonian accounts for various experimental 
data concerning the ferromagnetism and the transport 

• Curie temperature.i'@Ei) 
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• Resistivity (absolute value and magnetoresistance). 

• Spin and charge excitation spectra-SS 

• Scaling behaviors in charge and spin dynamics. 

The mechanism of MR in the single-crystal of these high T c a^mppjinds are understood from 
double-exchange alone, namely due to spin disorder scattering.BO'C^QEf 

There still exist controversies with some experiments and the DE model. An example is the loss 
of spectral weight near the FermiJevel observed by photoemission experiments,^- E$ the other is 
the optical conductivity spectrum© at uj ^ 1 e V which shows substantial loss of the Drude weight 
(integrated oscilator strength in the low frequency region) in contradiction with the specific heat 
behavior which does not show such a large mass enhancement. It should be noted that these 
experiments are highly sensitive to surface states, and there always remains an open question 
whether these experiments really jneasure the bulk states or are merely observing the surface 
states. Recently, Takenaka et alx3 reported that the optical spectrum of the "clean" surface 
prepared by cleaving shows larger Drude part compared to those prepared by surface filing, and 
is consistently large with respect to the specific heat estimate. Although the issue is still far from 
conclusive results, we should be very careful about such experimental details. 

Let us discuss the value of Jh in manganites estimated from these calculations. From exper- 
iments concerning spin and charge dynamics, data fitting gives Jh = 1-5 ~ 2eV. On the other 
hand, static and low frequency experiments such as T c and p(T) measurements are well reproduced 
by Jr = 3 ~ 5eV. Such a change in the value of Jh may be understood by the renormalization 
of interaction. While high energy experiments (such as &{uj)) directly measure the bare value 
of Jh, low frequency measurement is affected by the dressed quasiparticle involving the vertex 
renormalization. In this case, particularly, the effect of Coulomb interaction should be important. 
Since our DE Hamiltonian does not involve on-site Coulomb repulsion terms, elimination of the 
double occupancy of itinerant electrons are underestimated. This is compensated by increasing 
the value of Jh which also energetically prohibits the on-site double occupancy. Thus the vertex 
renormalization due to Coulomb repulsion may increase the value of Jh in the low-frequency scale. 

Discussion: Absence of polaronic behaviors in LSMO. Millis et al® discussed the 
discrepancies in DE model and experiments of LSMO with respect to (i) Curie temperature (T c ), 
(ii) temperature dependence of resistivity (p(T)), and (iii) absolute value of resistivity (p(T ~ T c j). 
They concluded that DE alone does not explain the MR of LSMO. 

What we have shown here is, however, that the accurate calculation of the DE model and the 
measurements of high quality crystal samples indeed give consistent results. Polaron effects, if 
any, should be small enough to be irrelevant. For example, Roder et alx$ discussed the decrease 
of T c due to the reduction of electron hopping caused by the polaronic effect. Thus, if polaronic 
interaction is included to the DE model hamiltonian, T c of the model becomes inconsistently small 
compared to the experimental value. 

Another point to be discussed is the experimental difference between LSMO and LCMO (see 
also §5.1| ). Since LCMO shows different behaviors in resistivity, e.g. temperature dependence of 
p(T) as well as the MR scaling relations, different physics must be taking place. If one assumes 
the polaronic arguments, crossover from metallic (LSMO) to insulating (LCMO) behavior with 
respect to the resistivity above T c could be understood in a following way. Insulating behavior 
might be due to small-polaron formation in systems with large electron-lattice coupling, while 
systems with smaller electron-lattice coupling might show large-polaron behaviors. 

However, this argument is in contradiction with experimental data. Pair distributionjunction 
(PDF) measurement of atomic displacements observed by inelastic neutron scattering^ shows 
that there exist large Jahn- Teller distortions in a dynamic way even in the metallic phase of 
LSMO above T c . As long as high T c compound LSMO is concerned, lattice distortion does not 
seem to affect transport properties as well as magnetic behaviors (e.g., T c ). Since large lattice 
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distortion is observed both in LSMO and LCMC^EjEj* polaronic discussion alone is not sufficient 
to explain the qualitative difference in p(T). 



5.3 CMR effects in low Curie temperature compounds 



What is the origin of the difference in resistivity behavior of LCMO compared to LSMO? 
As shown in §5.2, the DE model is sufficient to explain LSMO, and should be a good reference 
model to discuss LCMO. We discuss the half-metallic behavior in ferromagnetic states and the 
competition between ferromagnetic metals and charge ordered insulators. 



Theories of CMR based on the change in transfer integral. Several theories based 
on microscopic models have been proposed to explain the CMR phenomena in LCMO. Polaron 
effects of Jahn- Teller distortions are introduced to explain the metal-insulato^txansition from 
the point of view of large polaron to small polaron crossover by magnetism.KtS&J' The idea 
of AnHprson WaJizatioi] .due to spin disorder as well as diagonal charge disorder has also been 
discussedMBSMMM® 

Many of these proposals relate magnetism and transport through the change of the hopping 
matrix element 

t eS {T,H) occos(#/2), (79) 

where 9 is the angle between nearest neighbor spins, as discussed by Anderson and Hasegawa.S 1 
Namely, in general, these scenarios discuss the existence of the critical value of hopping t c . Change 
of the mean magnetic structure controls the value of t e g, and it is considered that some kind of 
transition occurs at t e g = t c . 

In a small polaron scenario, there exists a competition between the electron kinetic energy 
-Ekin oc t fT and the polaron energy E po \(oc A). It is characterized by the critical value for hopping 
t c where transition from small polaron state (small t e s region) and large polaron state (large t e ff 
region) occurs. Such a change in the behavior of carriers affect the conductivity through the self 
trapping mobility or the percolation of localized carriers. In the context of Anderson localization, 
critical value t c is determined by the relationship between the mobility edge and the Fermi level. In 
either cases, electron hopping t e s is considered to be controlled by the magnetic structure through 
the Anderson-Hasegawa's DE mechanism. In the presence of the competition between extended 
and localized states, metallic states are realized at the large hopping region t e s(T,H) > t c , and 
the system becomes insulating at small hopping region t ef[ (T, H) < t c . 

However, in order to explain the experiments for resistivity, these scenario need additional 
explanations to justify themselves. Let us point out two issues here: 

• Does the transition point t c really exist at the well-doped system x — 1/3 ~ 1/2? 

• Why does the "metal- insulator" transition occur only in the vicinity of T c , irrespective of 
carrier concentration and bandwidth? 

Localization phenomena in low carrier concentration is explained by the suppressed overlaps 
of wavefunctions. For example, in a low carrier small polaron system, a self-trapped polaron has a 
short confinement length scale compared to mean polaron-polaron distances. On the other hand, 
if we consider the region x ~ 1/3, every two sites out of total six neighbors in a cubic lattice 
is occupied by other polarons. In such a high density system, quantum mechanical overlaps of 
polaron wavefunctions have to be quite large. In other words, in order to assume a localized state, 
the gain of self-trapping potential energy has to be unrealistically large to compensate the loss of 
kinetic energy by such a localization with small length scale. For example, dynamic Jahn- Teller 
scenario requires that the polaron binding energy hfi larger than the electron kinetic energy to 
make insulating behaviors for the carrier doped case.EZP Thus, it is difficult to understand whether 
the localized state really exists in a realistic parametrization of the Hamiltonian. 
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Another point which is hard to explain is the fact that the metal-insulator transition always 
occurs in the vicinity of T c for various A-site substitutions in both ionic radius and average valence 
changes. In order to explain the metal-insulator transition scenarios require the "pinning" of the 
critical point at T ~ T c , i.e., irrespective of carrier concentration as well as A-site bandwidth and 
randomness, t c ~ t e s(T = T c ) always has to be satisfied. Note that t e g is determined by the short 
range spin correlation (Si ■ Sj) and is a smooth and continuous function of temperature without an 
anomaly at T c . It is also required to explain the complete absence of the metal-insulator transition 
in high-T c compounds such as LSMO at x ~ 1/3 within some realistic parametrizations. 



Magnetic inhomogeneities and the nanodomain TMR mechanism. In manganites, 
especially in the low-T c regions, magnetic inhomogeneities are experimentally observed. Here we 
discuss the importance of such inhomogeneities. 

j-llnconventional feature in the low-T c compounds is the presence of the central peak well below 
T C ,H23' ) which indicates the presence of the magnetic cluster and its diffusive dynamics. From the 
spin diffusion constant, the correlation length of the spin clusters are estimated to be £ ~ 10A. 
Neutron elastic ^scattering measurements also observed the ferromagnetic cluster with correlation 



length £ ~ 20A. 



104, 10 



Magnetic inhomogeneities also show up at the linewidth of the spin wave T q ~ g 2 ,0) which 
systematically increases as T c is suppressed by A-site substitution. It is speculated to be due to 
inhomogeneity effect through spin stiffness distribution.^'^ For lower T c compounds, there exist 
much prominent broadening of the spin-ypaae. ddspBisipn at the zone boundary!^ 

Optical conductivity measurements"'I^E^all2§tl2l' show that the A-site substitution causes 
changes in the spectrum at the peak structure around ~ 1.5eV and the infrared quasi-Drude 
(incoherent) structures. Formation of lattice polaron at lower T c compounds which causes uo ~ 
1.5eV peak is discussed. JLnhiirnjageneities in charge and magnetic structures observed by /iSR 



and X-ray measurements^' ti^liii' suggest that such polaronic cluster remain even at at low 



temperatures. Possible micrograin formation due to charge segregation as well as-nhase. separation 
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between ferromagnetic and antiferromagnetic domains has also been discussed.! 

Thus in most compounds inhomogeneous behaviors are observed experimentally. The tendency 
is that inhomogeneity is more prominent in the compounds with lower T c due to smaller (r (or 
larger a{rx) as well). 



Let us now focus on the case x ~ 1/3. As discussed in § |5.1| , decrease of T c is much larger 
than the estimate from the reduction of the bandwidth by (ta)- Inhomogeneity may be playing 
an important role to this behavior. Since wide bandwidth compound LSMO shows a ferromag- 
netic metal phase while narrow band compound PCMO is a charge-ordered insulator, there exist a 
competition as well as a bistability of these phases in the intermediate bandwidth region. Inhomo- 
geneous behaviors in LCMO suggests that coexistence of microscopic domains with ferromagnetic 
and charge-ordering correlations might happen. In the paramagnetic phase, neither correlations 
are long-ranged in a macroscopic sense. They should be short ranged and/or dynamic. 

In such nanodomain structures with microscopic phase separation, magnetic phase transition 
occurs when inter-domain correlations become long ranged. Although intra-domain correlations 
may begin at higher temperature determined by the DE mechanism, true long range order is 
controlled by the domain-domain interactions mediated by the junction structures. Intermediate 
region with charge ordering reduces the magnetic coupling between ferromagnetic nanodomains. 
Then T c should be substantially reduced. Such nanodomain structure creates two different energy 
scales for magnetisms, i.e. intra-domain and inter-domain interactions. 

It also shows up in two correlation lengths. One is the intra-domain correlation length, which 
becomes the domain size at low temperature regions. The other is the inter-domain correlation 
length which is the length scale to determine the nature of the magnetic phase transition. It has 
been reported that in the low T c compounds the magnetic correlation length does not diverge 
at T c but stays constant, e.g. ~ 20A for NSMOl^ This unconventional behavior is understood 
if the typical length scale of the nanodomain structure is ~ 20A. Due to the resolution of the 
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triple-axis experiment, it seems that the macroscopic correlation length which should diverge at 
T c was not detected. 

Now we discuss the relationship with CMR phenomena. In the presence of nanodomain struc- 
tures with ferromagnetic metals and (charge ordered) insulators, several mechanisms create mag- 
net oresistance. A possible mechanism is the percolation of metallic nanodomains. If the external 
magnetic field is applied, the system gains energy by increasing the volume fraction of the fer- 
romagnetic nanodomains. At the percolation threshold, there exists a metal-insulator transition. 
This scenario is, however, unlikely in the sense that metal-insulator transition occurs only at the 
percolation threshold, while the CMR phenomena is widely observed for various composition range 
away from some critical point. 

Another possibility tp_be discussed here is the nanodomain TMR phenomena, as is commonly 
the case for polycrystalslHf At T ~ T c , the intra-domain spin correlation is well developed so each 
nanodomain can behave as a half-metallic domains. Then, application of magnetic field controls 
the spin valve transports between nanodomains. This scenario is most likely in the sense that 
it naturally explains generic MR behavior at T ~ T c and does not assume any critical points. 
The conductivity is controlled by magnetism through the spin valve channels. Huge sensitivity to 
the external field is due to the fact that each metallic nanodomain already forms a ferromagnetic 
cluster. The idea is consistent with the phenomenologicaLexplanation of resistivity in LCMO 
by the Two-fluid model proposed by Jaime and Salomon.!^ They discuss the coexistence of a 
metallic conductivity path and an activation-type polaronic conductivity. 

Let us finally mention the origin of such nanodomain structures. Instability of electronic phase 
separation^ is a candidate for the initial driving force for such phenomena, which is stabilized 
to form droplet structures due to long range Coulomb interactions. Another possibility is the 
effect of static potential disorder due to A-site cation R 3+ -A 2+ distributions which causes charge 
inhomogeneities,tl3^ as well as self-trapping effect of lattice polarons. 



VI. SUMMARY AND CONCLUDING REMARKS 

In this article we discussed the thermodynamics of the DE model. We treated the model at 
finite temperature non-perturbatively with respect to the spin fluctuations. Various thermody- 
namic quantities including magnetic and transport properties are calculated. Modifications to the 
mean-field type treatment by Anderson- Hasegawa and de Gennes are made. 

These new results are found to be very important in comparing the experimental data for 
manganites with the DE model. As long as the high T c compounds (e.g. LSMO) in the metallic 
phase are concerned, the DE model accounts for various experimental properties including mag- 
netic transition and resistivity. Comparison with respect to low T c compounds (e.g. LCMO) are 
also discussed. 

There .still remains many open j^nest ions. One is the roles of other interactions such as lattice 
distortions' or orbital fluctuations.!^ HlJ It is experimentally clear that these are quite important 
and indeed make some long range orderings in the insulating phases. However, as we have shown 
in this article, such interactions does not show up in thermodynamic properties in the high T c 
metallic phases. It is quite interesting how and why such "screening" of interactions occur in the 
metallic phase. 

In other words, if we start from a realistic model for manganites with orbital and lattice degrees 
of freedom as well as Coulomb interactions etc., we somehow end up with the double-exchange 
model in the metallic phase as a consequence of neglecting high-energy excitations. Therefore, we 
should always regard the model as the renormalized model with parameters t = t e s and J = J e g 
in eq. (|3]). Or, in a strict sense, we should consider the action of the double-exchange model in 
eq. (0) with renormalized Green's function G c s and the coupling strength J e g. Such a detailed 
renormalization studies will help us to understand other intermediate phases. 

Another point of interest is to understand the nature of the inhomogeneous systems with 
nanodomain structures. Such study with respect to both microscopic and mesoscopic length scale 
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might help us to understand the generic features of the physics in strongly correlated oxides. 
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